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We present a general formula for tunneling conductance in ballistic ferromagnet/ferromagnetic 
insulator/superconductor junctions where the superconducting state has opposite spin pairing sym- 
metry. The formula can involve correctly a ferromagnetism has been induced by effective mass 
difference between up- and down-spin electrons. Then, this effective mass mismatch ferromagnet 
and standard Stoner ferromagnet have been employed in this paper. As an application of the 
formulation, we have studied the tunneling effect for junctions including spin-triplet p-wave super- 
conductor, where we choose a normal insulator for the insulating region whereas our formula can 
treat a ferromagnetic insulator. Then, we have been able to devote our attention to features of 
ferromagnetic metal. The conductace spectra show a clear difference between two ferromagnets 
depending upon the way of normalization of the conductance. Especially, a essential difference is 
seen in zero-bias conductance peaks reflecting characteristics of each ferromagnets. From obtained 
results, it will be suggested that the measurements of the tunneling conductance in the junction 
provide us a useful information about the mechanism of itinerant ferromagnetism in metals. 

PACS numbers: 72.25.-b, 74.25.F-, 74.20.Rp 



INTRODUCTION 



Andreev reflection^ AR), which occurs at the inter- 
face of the junctions involving superconductors, is one 
of the most important elemental processes in the trans- 
port through the superconducting junctions A the- 
ory of transport taking into account the AR was for- 
mulated by Blonder, Tinkham and Klapwijk referred 
to as BTK theory Q . The BTK theory enables us to 
probe the pairing state of superconductors. For exam- 
ple, for a junction consists of a normal metal and an 
unconventional superconductor, the quantum interfer- 
ence effect between the injected and the Andreev re- 
flected particles, which feel mutually the different sign 
of superconducting pair potential through the scattering 
event at the interface of the junction, forms the so-called 
zero-energy Andreev bound states(ZABS)Hi]. Indeed, 
the ZABS originated from the <i-wave symmetry of su- 
perconducting pair potential have been observed as the 
zero-bias conductance peaks(ZBCPs) in tunneling exper- 
iment of high Tq cuprate superconductors according 
to theoretical prediction of Tanaka and Kashiwaya(TK) 
formula3. The ZBCPs reflecting the ZABS arc the es- 
sential feature of the electrical conductions in normal 
metal / insulator / unconventional superconductor junction 
and provide us the information of the pairing symmetry 
of superconductor On the other hand, in ferro- 

magnct/insulator/conventional superconductor junction, 
the measurements for low energy transport via AR offer 
the opportunity to probe the magnetic property such as 
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the polarization of ferromagnetic materials [11|-[13[. The 
AR in this junction is suppressed by the exchange field 
in the ferromagnet layer. As a result, the conductance 
at low energy of the junction is suppressed respond- 
ing to the polarization of ferromagnet. The behavior 
of ZBCPs in fcrromagent /insulator/unconventional su- 
perconductor junction have been studied to understand 
the characteristic properties of unconventional supercon- 
ductors and to utilize the properties as applications of 
spintronics flil - flfjj . In these junctions, the ferromagnet 
has been described within the Stoner model based on a 
picture of free electrons. However, in some materials, the 
other descriptions of ferromagnetism are required. The 
ferromagnetism kinctically driven by a spin-dependent 
bandwidth asymmetry, or, equivalcntly, by an effective 
mass splitting between f- and .[-spin particles [i~7rl22| is 
an interesting model giving itinerant ferromagnet. 

Recently, Annunziata, ct al. analyzed 
the charge and spin transport in ferromag- 
net/insulator/superconductor(F/I/S) junctions with 
taking account above mentioned the spin-dependent 
bandwidth asymmetry ferromagnet (SBAF) making the 
effective masses have different values in ferromagnet [23l— 
[25| . They clarified that from the knowledge of the 
critical transmission angle the measurement of the 
effective mass difference among the particles could be 
possible [HI. In there, it has been also shown that the 
F/I/S junction is an effective probe to investigate the 
mechanism of ferromagnetism and the pairing symmetry 
of unconventional superconductor. Furthermore, it is 
suggested that the F/I/S junction can be useful as 
switching device using the spin current in the case of 
the symmetry of superconductor is conventional s-wave 
case[24| and be possible as a spin-filtering device[2Fj|. 
They have studied on F/I/S junctions with several types 
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of superconducting symmetries as the conventional s- 
wave, unconventional d a .2_ y 2-wave, and the time reversal 
symmetry broken d x 2_ y 2 + is or d x 2_ y 2 + d xy states. 
Moreover, although zero bias anomaly of differential 
resistance has been shown experimentally in Sr2Ru04-Pt 
point contact experiment prjl] . more recently, Kashiwaya 
et. al.[27| has shown the ZBCP which is direct evidence 
of the ZABS by tunneling spectroscopic experiment of 
Sr2Ru04-Au junction. There is much interest and im- 
portance to investigate problems on spin triplet p-wave 
symmetry nature because the p-wave pairing, especially, 
a chiral p x ± ip y state breaking time-reversal symmetry 
is one of the best candidates for bulk superconducting 
state of Sr 2 Ru0 4 HJIp. 

In this paper, a formulation of the tunnel- 
ing conductance for charge and spin currents in 
ferromagnet / ferromagnetic- insulator /superconductor 
(F/FI/S) junctions will be presented by taking the ef- 
fective mass difference leading the spin-band asymm etry 
between f- and |-spin particles in ferromagnet |23l - [25| 
into our previous theory Although the formulation 
can be used for singlet and S z = triplet superconduc- 
tors, we will study a chiral p-wave state. Our formula 
has general form being able to include the ferromagnetic 
insulator, however, a normal insulator surrounded by the 
ferromagnet and the superconductor is considered for 
the insulating layer to get pure characteristic features of 
a ferromagnetism in a ferromagnet. It is found that the 
normalized conductance spectra shows a clear difference 
between Stoner and spin-band asymmetry ferromagnets 
(STF and SBAF for abbreviation). The present results 
may be helpful in investigations of the mechanism of 
ferromagnet. 

This paper is organized as follows. In Sec. II we ex- 
plain a theoretical model and derive a formulation fol- 
lowing our previous method based on the BTK theory. 
The results for ferromagnet /insulator/ chiral p-wave su- 
perconductor junctions are presented in Sec. HI. Finally 
the results arc summarized in Sec. IV. 
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II. 



MODEL AND FORMULATION 



For the model of formulation, we consider a two- 
dimensional ballistic F/FI/S junction with semi-infinite 
electrodes shown in figure 1. A flat interface is assumed 
to be located at x=0, and the ferromagnetic insulator 
for up(down) spin is described by a potential Vfn\(x) 
= (Vb + (— )V ex )S(x), where S(x), Vq and V ex are the 5 
function, a nonmagnetic barrier amplitude and a mag- 
netic barrier amplitude, respectively. 

For the ferromagnetism in the F electrode, we adopt 
two kinds models of mechanisms shown in Fi^2] One 
of these is the standard Stoner model in which the fer- 
romagnetism is induced by the exchange potential lead- 
ing to the rigid energy shift between f-spin and J.-spin 
bands. The other is a spin bandwidth asymmetry model 
proposed by Hirsch[l8j], in which the bandwidth is tuned 
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FIG. 1: Schematic illustration of scattering processes of an 
injected electron with f-spin at the F/FI/S ballistic junction. 
Here, 9^, Oar, and 8 s are injection, Andreev reflection, and 
transmission angles, respectively. It is assumed that the nor- 
mal reflection at the interface is totally specular, then normal 
reflection angle is also given by 0^. For the case of 4--spm 
electron, it can be depicted by flipping ~\ by j, in the figure. 



relatively by the ratio of the effective masses between f- 
and 4--spin particles. Although in the following the free 
particle-like spectra of parabolic type is assumed as nor- 
mal electronic dispersion relation, we suppose to define 
the concept of bandwidth for the above description. It 
implies that there can be some relations between this 
description and some effective one-band tight binding 
model permitting the effective masses of carriers being 
proportional to the inverse of the width of the bands 
where the carriers get itinerancy. Hence, only giving dif- 
ferent values of the masses for and 4,-spin electrons 
yields a bandwidth asymmetry. 

The spatial dependence of the pair potential is taken 
as A(r) = AO(.t) for simplicity. In addition, we 
consider the S z = pairing states, where the ele- 
ments of pairpotential are given by A-j- f = A4.4 = 
and A^ = — Aj,^ for the singlet pairing state or 
A-)- 4 = A4.4 for triplet pairing state. Thus, the effec- 
tive Hamiltonian(Bogoliubov-dc Gnncs (BdG) equation) 
of the system can be reduced the decoupled equation for 
the eigen states (u F ^ sn (i)( r ), v F(s)i(t)(r)) T and is given 
by 
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Here E is the energy of the quasipaticle and Hq(t) is 
the single particle Hamiltonian for cr-spin where a = —a. 
In the Ferromagnet side, the single particle Hamilto- 
nian is given by H${r) = — ft 2 V 2 /2m a — pU ex — Efm 
where a =t, i, m a is the effective mass for cr-band par- 
ticle, p = +1(— 1) for | (Fj-spin, U ex is the exchange 
potential and Epm is the Fermi energy. The Hq in 
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FIG. 2: A sketch of dispersion relation between energy and 
wave number, and Fermi surface for STF (left side) and for 
SBAF (right side). It is assumed that free electron model 
with rigidly energy shift in STF and different effective masses 
for each spin in SBAF. 

the superconductor side is given by Hq(t) = H^(r) — 
— ft 2 V 2 /2ms — Eps where the m s and Eps are the ef- 
fective mass of the quasiparticle and the Fermi energy, 
respectively. To describe the Fermi surface difference, we 
assume Epm = Eps = Ep. 



In the following, we apply the quasiclassical approxi- 
mation where Ep ^> ( E, A(fe) ) and the fc-dependence 
of A(fc) is replaced by the angle 9s between the direc- 
tion of the trajectory of quasiparticles in the supercon- 
ductor and the interface normal. In the quasiclassical 
approximation, the wave vectors of and kpLQ{HLQ) 
are given by /c t(i) = y / (2m t(i )/h 2 )(E F + (-)U ex ) and 

kpLQ(HLQ) = k s = ^2m s Ep/h 2 , respectively, where 
ELQ (HLQ) indicates electronlike (holelike) quasiparti- 
cles. For example, we assume the injection of f-spin elec- 
trons from the ferromagnet at an angle (9-j- to the interface 
normal as shown in FigJTJ There are four possible scatter- 
ing trajectories exist; Andreev reflection with angle 9ab, 
as holes belonging to J,-spin band, normal reflection(NR), 
transmission to superconductor as ELQ, and transmis- 
sion as HLQ. These four processes are described in same 
way for .[-spin electrons with changing the scattering an- 
gle 9^ to Since the translational symmetry holds along 
the y— axis, the parallel momentum components of all 
trajectories are conserved k a sin 9 a =k s sin 9 S =kg sin 9s- 
The angles 9 a . 9 3 and 9s differ from each other except 
when U ex = and m^=mi=ms, which means retrore- 
flectiverly of AR broken by the exchange field and the 
effective masses difference. The BdG equations are re- 
duced to the effective one-dimentional equation due to 
the translational invariance along y-axis of the Hamilto- 
nian. Thus, the solutions of the BdG equations for cr-spin 
electron injections are described as 
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, A+ = A(f5 s ), A-=A(n-9s) 



(5) 



where the probability coefficients a s , b a , c a , and d a are for AR, NR, transmission ELQ and HLQ. These coefficients 
are calculated from the boundary conditions at x = 0, 



U{v)p a{s) (x = 0) = U(v)sa(a)(x = 0), 
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As explained in our previous paper, the reflection pro- 
cess depends upon the size relation of the Fermi surfaces 
between FM and SC. In the following, we will consider 
a situation where fcj, < k$ < fc-j-, and m^/ms = ms/mj. 
with m-f > ms > mj.. Following the BTK theory with 
taking care of the probability conservation of quasiparti- 
cle flow, 

I K | 2 +^ | a s | 2 | c a | 2 | d a | 2 = 1 (9) 

V f,<T V f,° V f,V 

the conductance G^^ for the cr-spin charge(spin) cur- 
rent through the system can be calculated by 



particles in ferromagnet and v s + ^ — hks/ms is that 
of the ELQ (HLQ) in superconductor. It is much worth 

C(S) 

to note that our conductance formula G s \ is different 



from former works [23l - l25l ]. On our way of formulating 
the conductance as an extension of previous formulation 
[1, [HI to the present situation, it has to be needed for 
correct treatment of mass mismatch in a same metal that 
the cofneient of AR a s should be given by the ratio of 
group velocities rather than that of wavenumbers as a 
consequence of conservation law of particle flow. Using 
the obtained AR (a s ) and NR ( b„) coefficients in the 
same way as the previous paper [14[ based on the TK 
formula^, the charge(spin) conductance for each spin j" 
and i. can be formulated by 



where Vf. a = hk a /m a is the group velocities of the cr-spin 
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where 0c = cos 1 y7 ^ 2 (x — 1 + 7 1 ^ 2 ) ( or 
sin -1 \J 7~ 1 / 2 (1 — x) ) is the critical angle of the 
AR measured in the superconductor side. Here, 
X = U ex /EF(0 < X < 1) and 7 = mf/mj. > 1. In the 
above, Gn.<j corresponds to the conductance when the 
superconductor is in the normal state. We calculate the 
normalized conductance defined by 

G c(s )(eV) = g? d °s cos0 s (P t G c s f + P^gf ) 
H'l d0 s cos 9 S {P^G NA + PiG N>l ) 

(14) 



where the polarization P a for cr-spin is expressed as 

p t = + ^ P . = Lx* 

7(1 + X) + 1-X' 1 7(1 + X) + 1-X' 

It is noted in general that the normalized conductance 
will be defined alternatively corresponding to the actual 
experiments. 

Above formulas (2.11), (2.12), and (2.14) can re- 
produce former formulas of tunneling conductance for 
junctions including triplet superconductor (TS). For 
m,f = mi, these eqations coincide to that of STF/I/TS 



5 



TABLE I: Numerical values of the magnetization M, the nor- 
malized exchange interaction x = U ex /Ep, and the mass ratio 
7 = mf/mj.. For the value M = 0.25, there are two cases, 
one is pure STF, \ = 0.25 and 7 = 1, and the other is pure 
SBAF, x = and 7 = 5/3. All other values of M are in the 
same way except the case of M = 0. 



magnetization M exchange Int.x mass mismatch 7 

1 

0.25 0.25 1 

0.25 5/3 

0.5 0.5 1 

0.5 3 

0.75 0.75 1 

0.75 7 

0.99 0.99 1 

0.99 200 



SBAF 

1 - 




100 200 
1 



FIG. 3: The magnetization M as a function of 7 in the pure 
SBAF, and inserted panel is M as a function of x in the pure 
STF. 



j unction [14j, and for m-f — m.4. and U ex = 0, the conduc- 
tance formula for N/I/TS junctional is reproduced. 



III. 



RESULTS 



At first, we notice about the growth of the magne- 
tization M for STF or SBAF. Using the polarization 
Pfj , M is given by M = - Pj,. For pure STF case 
(7 = 1), the magnetization is equal to the magnitude of 
exchange splitting M = \ (= U ex /Ep). For pure SBAF 
case(x = 0), the M is given by M = (7~l)/(7+l). Thus, 
the half metal state in SBAF case is unphysical situation 
because 7 = 00. Figure 3 shows the M in SBAF case 
as a function of 7. It can be seen that the growth rate 
of M becomes very gradual over 7 w 50. From this, one 
can expect the clear differences of transport properties 
depending on M between STF and SBA near the half 
metallic limit. Hereafter, we call "strong ferromagnetic 
regime" as a region under and near the half metallic limit. 

In the following sub sections, we apply our conductance 
formula to ferromagnet /insulator /triplet superconductor 
(F/I/TS) junction (F referred to as STF or SBAF) where 
V ex = 0. As the pairing potential, a triplet p-wave state is 
employed by choosing A n (6 s ) = A it (6» s ) = A exp(i9 s ), 
Aft (0s) = Aj4(#s) = for opposite spin pairing. And in 
addition, we choose some sets of parameters (x, 7) giving 
the same M = {0,0.25,0.5,0.75,0.99} shown in Table [Q 
so as to get clear characteristics of each ferromagnets. 



A. Distinction between STF and SBAF 

To investigate a consequence of the different mecha- 
nism of the magnetization, avoiding any effects of the 
normal barrier we consider the highly transparent junc- 
tion in the metallic limit (Zq = 0). In this case, the 
normalized total conductances G^{eV) show same trend 
that conductance values inside the energy gap eV < Aq 
are reduced when the value of magnetization M is in- 
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FIG. 4: Normalized conductance spectra for the charge cur- 
rent Gj-ieV) in the metallic limit, Z — in (a), and in (b), 
height of conductance at zero energy is plotted as a function of 
the magnetization M for STF (red line) and for SBAF (black 
line). 



creased for both STF and SBAF (Figgfa)). It indicates 
that the retro-reflectivity of the AR is broken due to the 
induced M . However, the M dependence of reduction 
for Gj,(eV < Aq) is different for each of them. The dif- 
ference can be seen more clearly in the M-dcpcndencc 
of conductance values at eV = 0, GS(0), in Figgl(b). It 
is found that the suppression of Gg(0) for SBAF case 
is weaker rather than that for STF case without weak 
magnetization regime, 0.0 < M ~ 0.2 and at half metallic 
limit, M = 1.0. To clear the reason of different M depen- 
dence of conductances for SBAF case and STF case, we 
show the critical angle of AR as a function of M in Figj5] 
The 6c for both SBAF and STF cases decreases with in- 
creasing M. It is found that the difference between angles 
is getting larger from M ~ 0.2 to ~ 0.9, and converges to 
zero at M = 1.0. For nearly half metallic limit M = 0.99, 
the 9c is almost suppressed in the STF case, while there 
still remains in the SBAF case. The critical angles for 
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STF and SBAF are C 

and B c = cos" ^ \ t+m j ■ \— °"* v vt+m- 
rcspectively. Then, it is clear that 6c in SBAF case is 



^l-(f^) 1/2 (= sin-y 
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FIG. 5: Critical angles of AR as functions of M for STF 
(black line) and for SBAF (red line). As described in main 
text, these angles are determined by the conservation con- 
dition for momenta parallel to the interface, fcj, sin 6± = 
fcssin^s. Such as in the present model, i.e., the ferromag- 
netism is given by a mismatch in kinetic energy the AR crit- 
ical angle is determined by 6c — ds for the case of 6^ — n/2 
because of satisfying the condition fcj, < ks. 



larger than that in STF case for same M except non- 
magnetic state, M = and half metal state, M = 1. 
Consequently, as shown in Fig [J] the Gj~{eV < A) in 
SBAF case is larger than that in STF case. 



B. Ferromagnetic feature on ZBCP 

It has been shown theoretically that the ZBCP in 
F/I/S junction would be useful for measuring the mag- 
netization of ferromagnet fbH flo| . In here, we study the 
validity of the ZBCP for the distinction of ferromag- 
ncts. Figure [6] shows the conductance Gj-(eV) for the 
junction in the tunneling limit Z = 5. The ZBCPs 
seen in both STF and SBAF cases are attributed to the 
anisotropy of the pair potential of p-wave superconduc- 
tor. For STF case, the previous results have been 
rcproduced(FigjBJa)). In contrast, there are some differ- 
ences for SBAF case. Especially, it is found that the con- 
ductance near eV = increases slightly with increasing 
M (Fig(B][b)). This opposite behavior can be seen more 
clearly in the M dependence of ZBCPiFigJBJc). With in- 
creasing M, in contrast to the monotonically decreasing 
behavior of STF case, the ZBCP in SBAF case increases 
up to a certain value of M in strong ferromagnetic regime 
and then, suddenly decreases toward the half metallic 
limit where the ZBCPs in both cases are suppressed per- 
fectly. Cause of this opposite behavior could be reduced 
to the definition of normalization way since the magni- 
tude of ZBCP being a constant value in non-normalized 
case depends on the conductance as the superconductor 
is in normal state. There is other definition of normal- 
ization by using the AR critical angle measured in the 
ferromagnet side [Til. [Fo|. However, in that case, the AR 
critical angle itself depends on and is controlled by the 
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FIG. 6: Normalized conductance spectra of the charge cur- 
rent G%(eV) in the tunneling limit, Z = 5 for STF/I/TS 
junction (a) and for SBAF/I/TS junction (b). And in (c), 
height of ZBCPs is plotted as a function of the magnetization 
M for STF (red line) and for SBAF (black line). 



magnitude of M , as a results, even the normalization 
depends on M . Accordingly, in order to avoid the influ- 
ence of M, we alternatively calculate an angle averaged 
conductance defined as in the following, 



Qs = <9s,t + Qs,i 

_ /;/ 2 2 dfl sCOS fl g p g Gg, e 



C/2 d0 S COS 6*5 



We show the calculated results of the angle averaged 
conductance Qs in Figl?] which, in both STF and SBAF 
cases, show same tendency to decrease as increasing 
M(FigJTJa)). Similarly, the ZBCP is decreasing function 
of M(FigJTJb)). It is also shown that the reduction ratio 
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FIG. 7: Angle averaged conductance spectra Qs as a func- 
tion of magnetization (a) and ZBCPs vs. magnetization 
strength (b) in the tunneling limit Z = 5. Here, the con- 
ductance for SBAF case is indicated as solid line and for STF 
case is as dotted line. 
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FIG. 8: Magnetization dependence of angle averaged normal 
conductance Qn, Qiv.ti an d Qn,i for SBAF case (solid line) 
and STF case (dotted line) in the tunneling limit for the case 
of the superconductor is in normal state, (a) for Z = 0, (b) 
for Z = 5. 



differs in each of both cases as same as that in metallic 
limit. Thus, the opposite behavior seen in normalized 
conductance would reduce to the conductance in normal 
state. Therefore, it is noticed that the conductance of 
the junction for the superconductor being in normal state 
play an important role on our attention for two different 
ferromagnetisms . 

In order to clarify the difference between STF case 
and SBAF case more, we calculate the conductance in 
ferromagent/normal metal (F/I/N) junction for both 
in metallic and in tunneling limits. The angle aver- 
aged conductance in F/I/N junction Qn = ^ CT Qn,g 
is defined in similar way to that in F/I/S junction 
replacing Gg by G N . The calculated results of 
Qn.c for both Z{) = and Zq = 5 are shown in 
FiglU The angle resolved conductance G N a for cr-spin 

is rewritten by G Na = 4 cos 6 s X a I ((cos 9 s + X a ) 2 + 



Z% a ) where A 



v/cos 2 S + p X in STF/I/N and 



7" 



-p/2 



y/cos 2 6 S + {Y /2 - 1) in SBAF/I/N junc- 
tions. Here, we mention properties of M-dcpcndcncc of 
G N through \ or 7 in advance of descriptions about 



Qn- In STF/I/N junction, the G N * increases following 
growth of the magnetization, i.e., with increasing x since 
the gain of Fermi energy due to the band shift is larger 
than the Fermi surface effect [l4[ acting as an effective 
barrier between STF and normal metal, under the con- 
servation of the momentum along y— direction. On the 
other hand, because there is no Fermi energy gain from 
spread of the band width due to the effective mass mis- 
match in SBAF and the influence of the effective barrier 
arising from the Fermi surface effect becomes stronger 
with the increase of 7, the G N t in SBAF/I/N junction 
decreases with increasing 7 and become zero in the limit 
of 7 — > 00. G N , for both STF and SBAF cases decreases 
with increasing the magnetization caused by x or 7. 

In the metallic limit Zq = (FigEJa)), it is found that 
the Qjv,t m STF/I/N junction increases with increasing 
AI in contrast to Qn,]. decreasing toward zero in half 
metal state. In this case, M is given directly as M = x- 
Thus, the total conductance Qn = Qn.-\ + Qn.i is re- 
duced slightly by the Fermi surface effect with increasing 
M up to ~ 0.7. In SBAF/I/N junction, we can see simi- 
lar behavior in Qn^U) ■ The increase of Q jv,t is owing to 
P T which is an increasing function of 7. However, near 
the half metallic limit, Qn,^ reduces rapidly reflecting the 
behavior of Gjv T which is a decreasing function of M to- 
ward zero at M = 1(7 = 00) as mentioned above. Thus, 
as shown in FigJH the G^(eV) in SBAF/I/S junction 
decreases slowly with increasing M with comparing to 
that in STF/I/S junction. The difference between STF 
and SBAF becomes more clearly in the tunneling limit 
Z = 5(FigH{b)). With increasing M, Q N , a in STF case 
varies in rapidly rather than that in SBAF case. This is a 
difference of a barrier effect felt by particles with cr-spin 
in each cases. The barrier potential simply becomes rel- 
atively lower for particles with f-spin and higher for par- 
ticles with J,-spin in the STF case due to the rigid Fermi 
energy shift. However, the particles in SBAF directly feel 
the barrier potential because there is no shift of the Fermi 
energy. Thus, in SBAF/I/N junction, the increase of the 
magnitude of Qn,^ due to P T is suppressed by the Fermi 
surface effect and barrier potential and then, Qjy,t is get- 
ting lower with increasing M in contrast to the STF case. 
Therefore, the Qn in SBAF/I/N junction shows the op- 
posite behavior of that in STF/I/N junction. As a result, 
the normalized conductance Gf (eV) in SBAF/I/S junc- 
tion increases due to the reduction of the Q n depending 
on M (Figj6pb)-(c)). Indeed, as shown in FiglJl the angle 
averaged conductance Qs(eV)s for both STF and SBAF 
case show same trend on varying M. Thus, it can be 
conclude that the measurement of Qn will be also useful 
to identify the STF and SBAF. However, we emphasize 
that the measurement of ZBCP originated from ZABS 
is more powerful probe to investigate ferromagnet than 
that of Qn- Because, two Qn$ seemingly show drasti- 
cally different behavior depending on M for enough large 
Zo (FigHJb)), by carefully looking of the figure, differ- 
ences of each values of Qns are not so large for same M 
except strong ferromagnetic regime. Therefore, it seems 
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that an experimental distinction will become more diffi- 
cult on measurement of Q^. The ZBCP is getting more 
clear for larger Zq, then which can be expected to play a 
role of good manifestation of the difference of STF and 
SBAF. 



IV. 



SUMMARY 



In summary, we have derived a formula of the 
tunneling conductance in ferromagnet/fcrromagnctic- 
insulator/superconductor with antiparallel spin pair- 
ing junction by extending our previous theory for 
standard Stoner ferromagnet (STF) so as to in- 
clude spin-band asymmetry ferromagnet (SBAF) orig- 
inated from effective mass mismatch between parti- 
cles with opposite spins. Applying the formulation to 
ferromaget/insulator/p-wave superconductor junctions, 
differences between pure STF and pure SBAF have been 
investigated intensively. We found that, with growing 
the magnetization, the difference becomes clear in tun- 
neling conductance. The clarity of difference between 
STF and SBAF depends on the way of normalization of 
conductance and comes out more clearly in ZBCP near 



half- metallic limit. The obtained results suggest that the 
measurement of ZBCP may be useful for discriminating 
mechanism of ferromagnctism. 

Although our formulation includes the ferromagnetic 
insulator, we have studied only the normal insulating 
barrier case in this paper. The spin-filtering effect have 
been expected in the ferromagnetic insulator [f3| or in 
ferromagnet given by the effective mass mismatch[25|. 
Then, as an interesting future problem we will study ex- 
tensively the spin-filtering effect in junctions of including 
both ferromagnetic insulator and mass mismatch ferro- 
magnet connected to superconductors of s-, d-wave and 
broken time reversal symmetry pairing states. Moreover, 
it will be an important issue that the proximity effect is 
taken into account to the present formulation by carrying 
out the self-consistent calculation of the pairing poten- 
tial in order to analyze the actual experiments. Indeed, 
the ZBCP have been observed in tunn eling experiment 
of F/I/ii-wave superconductor iunction[35l|. And also, 
ZBCP in Si - 2Ru04 junction has been observed 27|, then, 
tunneling spectroscopy of F/I/Sr2Ru04 junction seems 
to be realized in near future. Our conductance formula 
can apply to such situations easily and get comparable 
results to experimental one. 



[1] A. F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964). 

[Sov. Phys. JETP 19, 1228 (1964)]. 
[2] G. F. Blonder, M. Tinkham, and T. m. Klapwijk, Phys. 

Rev. B 25, 4515 (1982). 
[3] C. R. Hu, Phy. Rev. Lett. 72, 1526 (1994). 
[4] Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, 3451 

(1995). 

[5] S. Kashiwaya, Y. Tanaka, M. Koyanagi, H. Takashima, 

and K. Kajimura, Phys. Rev. B 51, 1350 (1995). 
[6] S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63, 1596 

(2000), Phys. Rev. B 56, 7847 (1997). 
[7] G. Deutscher, Rev. Mod. Phys. B 238, 109 (2005). 
[8] M. Yamashiro, Y. Tanaka, and S. Kashiwaya, Phys. Rev. 

B 56, 7847 (1997). 
[9] C. Honerkamp and M. Sigrist, J. Low. Temp. Phys. Ill, 

895 (1998). 

[10] T. Lofwander, V. S. Shumeiko, and G. Wendin, Super- 

cond. Sci. Technol. 14, R53 (1997). 
[11] M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. 

Lett. 74, 1657 (1995). 
[12] R. J. Soulen Jr., J. M. Byers, M. S. Osofsky, B. Nadgorny, 

T. Ambrose, S. F. Cheng, P. R. Broussard, C. T. Tanaka, 

N. Nowak, J. S. Moodera, A. Barry, and J. M. D. Coey, 

Science 282, 282 (1998). 
[13] S. K. M. Upadhyay, A. Palanisami, R. N. Louie, and R. 

A. Buhrman, Phys. Rev. Lett. 81, 3247 (1998). 
[14] N. Yoshida, Y. Tanaka, J. Inoue, and S. Kashiwaya, 

J. Phys. Soc. Jpn. 68, 1071 (1999); S. Kashiwaya, Y. 

Tanaka, N. Yoshida, and M. R. Beasley, Phys. Rev. B 

60, 3572 (1999); 

[15] I. Zutic and O. T. Vails, Phys. Rev. B 60 (1999) 6320; 

61, 1555 (2000); J.-X. Zhu, B. Friedman, and C. S. Ting, 
Phys. Rev. B 59, 9558 (1999); J.-X. Zhu and C. S. Ting, 



Phys. Rev. B 61, 1456 (2000); Z. C. Dong, D. Y. Xing, 
Z. D. Wang, Z. Zheng, and J. Dong, Phys. Rev. B 61, 
144520 (2001); N. Stefanakis, Phys. Rev! B 64, 224502 
(2001); J. Phys. Condens.Matter 13, 3643 (2001); 

[16] T. Hirai, Y. Tanaka, N. Yoshida, Y. Asano, J. Inoue, 
and S. Kashiwaya, Phys. Rev. B 67, 174501 (2003); I. 
Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 
323 (2004); J. Linder and A. Sudb0, Phys. Rev. B 75, 
134509 (2007); P. H. Barsic and O. T. Vails, Phys. Rev. 
B 79, 014502 (2009). 

[17] C. Zener, Phys. Rev. 82, 43 (1951); P. W. Anderson 
and H. Hasegawa, Phys. Rev. 100, 675 (1955); P. G. 
de Gennes, Phys. Rev. 118, 141 (1960). 

[18] J. E. Hirsch, Phys. Rev. B 40, 2354 (1989); 40, 9061 
(1989); 43, 705 (1991); 59, 6256 (1999); 62, 14131 
(2000); Physica C 341-348, 211 (2000). 

[19] D. K. Campbell, J. T. Gammel, and E. Y. Loh, Phys. 
Rev. B 38, 12043 (1988); 42, 475 (1990); S. Kivelson, 
W.-P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 
Lett. 58, 1899 (1987). 

[20] Y. Okimoto, T. Katsufuji, T. Ishikawa, A. Urushibara, T. 
Arima, and Y. Tokura, Phys. Rev. Lett. 75, 109 (1995); 
Y. Okimoto, T. Katsufuji, T. Ishikawa, T. Arima, and 
Y. Tokura, Phys. Rev. B 55, 4206 (1997); S. Broderick, 
B. Ruzicka, L. Degiorgi, H. R. Ott, J. L. Sarrao, and Z. 
Fisk, Phys. Rev. B 65, 121102 (2002). 

[21] M. Higashiguchi, K. Shimada, K. Nishiura, X. Cui, H. 
Namatame, and Masaki Taniguchi, Phys. Rev. B 72, 
214438 (2005). 

[22] A. McCollam, S. R. Julian, P. M. C. Rourke, D. Aoki, 
and J. Flouquet, Phys. Rev. Lett. 94, 186401 (2005); I. 
Sheikin, A. Groger, S. Raymond, D. Jaccard, D. Aoki, 
H. Harima, and J. Flouquet, Phys. Rev. B 67, 094420 



9 



(2003). 

[23] G. Annunziata, M. Cuoco, C. Noce, A. Romano, and P. 

Gentile, Phys. Rev. B 80, 012503 (2009). 
[24] G. Annunziata, M. Cuoco, P. Gentile, A. Romano, and 

C. Noce, Phys. Rev. B 83, 094507 (2011). 
[25] G. Annunziata, M. Cuoco, P. Gentile, A. Romano, and 

C. Noce, Supercond. Sci. Technol. 24, 024021 (2011). 
[26] F. Laube, G. Goll, H. v. Lohneysen, F. Fogelstrom, and 

F. Lichtenberg, Phys. Rev. Lett. 84, 1595 (2000). 
[27] S. Kashiwaya, H. Kashiwaya, H. Kambara, T. Furuta, 

H. Yaguchi, Y. Tanaka, and Y. Maeno, Phys. Rev. Lett. 

107, 077003 (2011). 
[28] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. 

Fujita, J.G. Bednorz, and F. Lichtenberg, Nature 372, 

532 (1994). 

[29] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. 



Mao, Y. Mori, and Y. Maeno, Nature 396, 658 (1998). 
[30] G.M. Luke, Y. Fudamoto, K.M. Kojima, M.I. Larkin, J. 

Merrin, B. Nachumi, Y.J. Uemura, Y. Maeno, Z.Q. Mao, 

Y. Mori, H. Nakamura, and M. Sigrist, Nature (London) 

394, 558 (1998). 
[31] A.P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 

(2003). 

[32] K.D. Nelson, Z.Q. Mao, Y. Maeno, and Y. Liu, Science 

306, 1151 (2004). 
[33] M. Sigrist, Prog. Theor. Phys. Suppl. 160, 1 (2005) 
[34] Y. Maeno, S. Kittaka, T. Nomura, S. Yonezawal, and K. 

Ishidal, J. Phys. Soc. Jpn. 81, 011009 (2012). 
[35] A. Sawa, S. Kashiwaya, H. Obara, H. Yamasaki, M. Koy- 

anagi, N. Yoshida, and Y. Tanaka, Physica C 339, 287 

(2000) 



